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Matter gets organized at several levels of structural rearrangements. At a mesoscopic level one 
can distinguish between two types of rearrangements, conforming to different close-packing or den- 
sification conditions, appearing during different evolution stages. The cluster formations appear 
to be temperature- and space-dimension dependent. They suffer a type of Verhulst-like saturation 
(frustration) when one couples the growing (instability) and mechanical stress relaxation modes 
together. They manifest a chaotic behavior both in space and time domains. We pretend to offer a 
comprehensive and realistic picture of a material or mega-cluster formation in d dimension. 

I. INRODUCTION 

Matter organisations at a mesoscopic (molecular-cluster) level typically manifest a multitude of 
microstructural rearrangements. Cluster-cluster aggregations of proteins and/or colloids, phase sepa- 
rations, flocculation-coagulation phase transformations, sol-gel systems, (wet) sand or rice piles, etc., 
are manifestations of loosely-packed rearrangements, typically occurring under moderate or high tem- 
perature conditions. In contrast ripened polycrystals, sintered powders, soap froths and bubbles, and 
other cellular systems, constitute a type of rearrangement that usually emerges in a (relatively) low 
temperature limit and under certain ("field dependent") matter close-packing constraints. Beside such 
an agglomeration, fracture, desaggregation, desorption, dissolution, and alike, can be thought to be 
the "inverse process", finding its place in the opposite part of the relevant phase diagram ]]. 

In all of them spatial as well as temporal signatures of chaotic behavior, due to matter reorganisa- 
tions, can be detected: They are temperature and space-dimension dependent. In particular, one can 
show rigorously that in the limit of the spatial dimension going to infinity loosely-packed agglomera- 
tions become non-chaotic by suppressing totally their instability growing mode since it is related to the 



noncquilibrium agglomerate's entropy Ql, while their closed-packed counterparts are not. When the 
growing mode is coupled to a mechanical stress relaxation mode as a power law via some phenomeno- 
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logical relation of Hall-Petch-Griffith (H-P-G) |3 type (an Onsager-type conjecture |^] of the present 
study), certain marks of Bethe-lattice frustration, related to a spatial overcrowding of the Cayley-tree 
branches, appear in the (mean-field) approach - a kind of frustration qualitatively of very similar type 
than that observed in Verhulst-type systems in an adequate time and parametric zone. 

In the specific cases discussed in this review, however, by increasing the space dimension, d, we 
automatically induce some increase of the possible number of degrees of freedom in the system. Thus, 
when taking into account the coupling of the late-time growing and relaxation modes, say, in a fairly 
synchronized viz power-law way of H-P-G type, one obtains that even though the material's relaxation 
goes slower than in the case when such a coupling is proposed in an unsynchronized (Debye-relaxation 
involoving, i.e. rapid) way, one is, however, able to establish or restore an apparent dynamic mi- 
crostructural order within the system the noncquilibrium (chaotic) measures of which are proposed 
below. There is, unfortunately, no way of establishing such an order when the coupling fails the 
power-law type synchronization requirement 

Thus, when an ample space amongst the clusters is recovered by the system at its mature grow- 
ing stage, we consider that the system successfully tries to avoid a chaotic matter organisation in 
space. Note that temperature may markedly help in surmounting the activation-energy barrier of 
the agglomeration, especially when it is raised appropriately, whence when not "damaging' a possibly 
smooth evolution of the system. Full success is, however, guaranteed when the limit of d ^ oo is 
reached. If there is no chance for recovering the ample space, the late-time growing stage is realized in 
a moderately chaotic way. The mechanical stress relaxation, in turn, enters a readily chaotic regime, 
since the (noncquilibrium) entropy of the system diverges to plus infinity. The overall scenario resem- 
bles, in general, a formation of large (fractal) colloid aggregates that typically occurs with and without 
temperature and/or space-dimension dependent gravity factor domination, like as if imposing some 
limits to gelation of colloids ^. 

The paper is arranged as follows. In Section 2, we define both the closely-packed as well as loosely- 
packed agglomerations, calling the latter the aggregation throughout. In Section 3, we list some 
qualitative signatures of chaos in matter-agglomerating systems, and refer briefly to different definitions 



as well as meanings of chaos. In Section 4, we present quantitative measures of chaos signatures in 
systems of interest, whereas in Section 5 we unveil number-theoretic measures, featuring a chaotic 
spatio-temporal behavior of them. In Section 6, according to some suggestions given in on which 
much of our report is based, in order to see which agglomerations behave orderly or non-chaotically, 
we explore the limit of d oo, and arrive at a certain interesting (perhaps, surprising) conclusion, 
favoring aggregation of matter, or some structural loosely-packed, and typically high-temperature, 
matter rearrangements - in contrast to those emerging under close-packing low-temperature conditions. 
We close the paper by offering a concluding address in Section 7. 



II. AGGLOMERATION VS. AGGREGATION OF MATTER - A MODEL DESCRIPTION 

Following 7] , throughout the present study, we wish to distinguish between the notions of agglomer- 
ation and aggregation of matter. By the former we mean an assembly of grains or molecular clusters, 
kept together by relatively strong forces {e.g. ionic), so that there is no easy possibility of taking the 
clusters apart, or destroying them. For the latter, because of the appearance of weak bonds between 
clusters, such as Van der Waals or hydrogen types, the possibility of cluster separation becomes an 
observed tendency of the matter rearrangement due to their weak bonds. For a schematic explanation 
of the difference between both matter arrangements, see Fig. 1. 



A. Basic system of equations describing model matter agglomeration 



As in previous work 



we begin with a local continuity equation 



(1) 



preferrentially supplemented by the corresponding initial (of delta-Dirac type as a first attempt) and 
boundary (typically, of absorbing type) conditions (IBCs). 

In a few subsequent studies a thermodynamic-kinetic description ^ ^, ^ of model complex matter 
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agglomeration has been worked out. For the current |32J in the space of cluster volumes 



Jiv, t) = - Biv) — ^\ fiv, t) - D{v)^f{v, t), 



(2) 
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FIG. 1: Typical cluster-merging (three-grain) scenario for closely-packed (left, denoted by A) and loosely- 
packed (right, denoted by B) agglomerations. Two consecutive time steps ti and t2 are shown. The former 
usually goes by a scenario with the preservation of the total agglomerate's volume (though in a more irregular 
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way, when its logarithmic speed is measured, cf. Sections 3-6), whereas the latter does not H. In the former, 
the clusters do not perform a translational motion but their boundaries may fluctuate in time and space, even 
though they are quite strongly confined by their neighborhood. In the latter, an almost opposite situation in 
the time-and-space domain is typically observed. Some void is left behind a loosely packed system 



has been used [9|, where f{v,t) is the distribution of clusters of volume v: this means, that f{v,t)dv 
is the (relative) number of clusters with size in the infinitesimal volume interval [v,v + dv]; t is the 
time; $ represents the physical potential, equivalent to the free energy of the system (see joj for an 
explanation of the term). It is assumed to be one of the most relevant drivers of the agglomeration 
process at the mesoscopic level, assuring its nonequilibrium character. In fact, the current ©J in 
the form presented above, comes from a rigorous derivation, starting from the Gibbs equation for the 
entropy production |ic| . It has quite strong foundations anchored in nonequilibrium thermodynamics 



It is worth mentioning that upon inserting Eq. Q into the continuity equation, Eq. JQ), one gets 
a second-order partial differential equation of the Fokker-Planck-Kolmogorov (F-P-K) type [l^ . 



For the dynamics of such a system some routes to chaotic behavior have been sketched elsewhere 
by considering the (in)stability of Markov semigroups in ,12iJ. The mobility B{v) is also defined in the 
configurational w-space and reads 



iviaj 

a 



where _D is a diffusion reference constant. Realize that the principal role of D is to scale the time 
variable; /cb represents the Boltzmann constant. Note that the mobility B{v) is related to the On- 
sager coefficient, L{v), that appears in the derivation of the matter flux equation Q under a set of 
assumptions, mostly based on the locality of the Kramers-type process [lol etc. L{v), and hence 
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B{v), could be measured by comparing the current and the thermodynamic force 



The quantity 



D{v) — Dv" is to be inferred from the Green-Kubo (G-K) formula Il4j, so that there is some quite 



strong suggestion for deriving B(v) both, experimentally as well as 



theoretically 



|. For instance, it is proved 



There is a debate about a possible violation of the G-K formula 
that for a gas of charged particles subjected to an external electric field, the mean mobility of a charged 
particle, based on the G-K formula, is reliably well estimated for suitably small values of the external 
field. Moreover, at a microscopic scale one observes a nonlinear (or chaotic) behavior of the particles, 
which is, unfortunately, not reflected by the macroscopic (mobilility) measure. In our case, we assume 
algebraic correlations in u~space, for a phenomenological formula. The assumption seems to be as 
natural as possible: D(v) oc v"; that means that both the diffusivity D{v) and the mobility B{v) 
are proportional to the cluster hypersurface, s'''' R'^^^ . It should be underlined that it is, in our 
opinion, the common physical case in clustering phenomena, and is working properly at the mesoscopic 
level considered in our approach. 

Notice right here that exactly the same assumption has been used to model in a F-P-K way the 



formation of surface nanostructure arrays 



Therein, an experimentally-observed passage between 
direct curvature-dependent ripening of matter nano- islands (our densely-packed agglomeration), and 
inverse ripening, with an elastic-field caused contraction of growing quantum dots '17'| (our sparsely- 
distributed agglomeration of matter), has been presented. 

There are, however, matter agglomerations, for a given T, that do not conform usually to 



D{v) oc B{v) ~ s^'^\ (4) 

To them belong both some physical- metallurgical transformations [iJ l. such as martensitic, and pre- 
sumably, also certain phase orderings of non-diffusive kind, emerging in model biosystems, such as 
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those occurring in lipid biomembranes |19l] . Other than algebraic types of correlations in the hyper- 
space can likely be expected for these. If a power law of the type given by Eq. Q can be kept for 



further modeling, some additional correlations in time must complete a more comprehensive correla- 



tional proposal, cf. |2C|. Other types of correlations in the hyperspace, even if they allow to get a 
general solution to the problem, may not accommodate the boundary conditions 9], so that one would 
expect either to be left with an unsolved specific problem or to encounter anomalous or irregular be- 



havior of the agglomerating system 



boundary boundary conditions 



2] 



2CM . In such a case another type of finite, instead of infinite 

n 

, can sometimes give a remedy for the problem 21]. Here, under 
the term infinite boundary conditions we typically understand the boundary conditions of absorbing 
(Dirchlet) type 



f{v = 0,t) = f{v^Vc,usut)^Q 



(5) 



in which the single cluster volume is taken at infinity, Vdust = 00, whereas in case of the finite 
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boundary conditions it assumes a finite value, < Vdust < 00, cf. 2l], and a discussion therein. 
Although the latter unquestionably seems more physical the former is more frequently used to reveal 



the evolutions in matter-agglomerating systems 



this resembles to some extent a situation in 



statistical-thermodynamical systems undergoing an equilibrium phase transition: As such they are 
typically considered in the so-called thermodynamic limit (here, with a number of subunits going to 
c») under the mentioned agglomeration-oriented, e.g. condensation conditions, and the analogy would 
presumably extend over the examined nonequilibrium evolutions too ^ ^ . 



B. Thermodynamic potentials driving matter agglomerations 



In previous work 9( the analytic form of a (so-called) compaction potential was obtained, i.e. 



(6) 



where $0, Ro - constants, and R stands for some cluster radius. Because 



V = V 



^'''>r^R'^, d=l,2,3,... 



(7) 



one gets also $(w) cx ln{v/vo), where Vo is a constant. 

The logarithmic potential assures the emergence ofjathcr compact and curvature-involving struc- 
tures, whence the name of "compaction potential" _22]. It should be noted that (f> is an entropic 



potential 9J. Thus, it can be a cause of some desaggregation, or matter-influenced impingement 
effects, occurring within the overall aggregation space. 

In a previous study on the phase transformation kinetics for loosely packed " diffusive" agglomerates 
we have written the matter flux of a purely diffusive nature prescribed in configurational space as 
follows 

J(v,t) ^ ~Div)^fiv,t). (8) 

(The diffusion function D{v) = Dv"" is proportional to the cluster (grain) surface.) 

Both closely-packed and loosely-packed agglomerations follow from the general form 0. Indeed, 
the loosely-packed case is obtained when the first (drift) term in r.h.s. of Q can be neglected. 



Formally, B{v) — > when T — > oo. From the physical point of view, it corresponds to 



sufficiently high 



above which 



temperatures T > Tpass > 0, where Tpass can be treated as a cross-over temperature 
the agglomeration takes place exclusively by yielding loosely-packed microstructures. However, the 
drift term in r.h.s. of Q depends both on v and T. It tends uniformly, which means independently of 
u, to zero at the high temperature limit if 

^ 5$(t^ ~ C = const. (9) 
ov 

Then for a given system, temperature Tpass does not depend on v and looks consistently defined. 

Some additional argumentation can be provided that such a constant (limit) C exists and is well- 
defined. Namely, when applying both ||2Jl and (T)) one sees with sufficient accuracy that 

This means that C is essentially determined by a product involving two contributions: a certain 
curvature-like term, k — 1/ \ AR |, and some dimensionless energetic argument, eg = | A$(w) j/fesT. 
The above claimed high-temperature limit, with the cross-over temperature Tpass as a reference tem- 
perature characteristic of a system of interest, would naturally demand < eg << 1 while, because of 
approaching the mature growing stage any change in the cluster radius must be small, <| AR |<< 1, 
and therefore, its inverse would tend to some big value, i.e. k >> 1. Thus, C will take on a finite 
value. It is believed that for certain agglomerations under readily high-temperature conditions it will 



In 



a condition of setting the current equal to zero, J = 0, has been chosen to balance diffusive 



and non-diffusive terms in the F-P-K type description, cf . [lOj , aiming at getting a proper behavior of 



eventually acquire a small valuels^, that means, < C << 1 naturally holds, ft is a case when the 
potential 

cx ~ (11) 

the metastable nanostructure [17 1 arrays. We are of the opinion that such a proposal is legitimate in the 
relatively low-temperature domain. When the temperature is raised, but agglomeration is still allowed 
to occur, the proposal may fail. Thus, the above is a possible solution for the high-temperature limit. 
A type of localization of the Gaussian distribution, characteristic of the inverse ripening (a metastable 
state of the nanostructure evolution) can also be obtained within the present modeling, cf. |^ . This is 
the case of Eq. ||SJ) when in a (readily) mature growing stage, since the single volume v of the cluster 
does not change much. As a matter of fact, there is no small-cluster population available for merging 
(Fig. 1), i.e. D(v) — > const., which nearly corresponds to the high-temperature criteria of Eq. @, or 
equivalently Eq. H1U|I . In so doing, Eq. (|SJ| represents the 1st Fick law in its standard form. Upon 
inserting it into Eq. ^ one immediately arrives at the 2nd Fick law (in the configurational space) 
with its standard Gaussian solution, the metastable case being emphasized in 1^ . 

The above potential form (|ll|l . designed for loosely-packed agglomeration, seems legitimate here: 
Note that the 'force' Fc-c oc d^{v)/dv behaves like 



- ^ - i' (12) 

because v ~ W^. Thus, Fc-c acts as the inverse of the area of the cluster hypersurface, s^''^ which 
implies that the smaller the area is, the bigger the force acting on the cluster can be, this way impeding 
the formation of new clusters, which would contribute to an aggregate's density increase. Qualitatively, 
a similar dependence is found for the closely-packed matter agglomeration: from ^ one gets, as above, 
for the 'force' 

^^-^ - (13) 
Here, Fc-c acts as the inverse of the hypervolume of the cluster, v^'^^ := R^, which makes a clear 



difference between ciosely-paclced and ioosely-paclced agglomerations, presumably leading to a certain 
relaxation of the surface tension conditions for loosely-packed clusters-containing systems 2^ . 
Referring further to Q and using the similarity relation, Eq. Q, one gets 

<i>(i?)(x$(i?„)^, (14) 
Ho 

where 



HRo) = ^Ro, T > Tpass, (15) 



and consequently, ^{R) oc R. Moreover, 

D„ = D(l-a). (16) 
Ro can now be specified to be the initial cluster radius. Note that Da is a d-dependent quantity. 

C. Cluster volume fluctuations as reliable characteristics of matter agglomeration 

Aggregations and agglomerations emerge in a fluctuating changing medium. Therefore, any reason- 
able quantitative attempt on resolving the fluctuation impact on their speed is worth examining here. 
In what follows, let us propose an evaluation of the reduced variance 

a^it) = <-'i^)>-<-y) > , < > _ 1, (17) 

as a direct measure of the cluster volume fluctuations. 

The notation used in Eq. H17|l refers to the statistical moments 

/•oc 

<i;«(t)>=/ v''f{v,t)dv n = 0,l,2,... (18) 

"'0 

of the stochastic process, where the matter agglomeration is usually described by the local continuity 
equation, Eq. |Q. 

The explicit solutions, f{v,t)s, have been presented elsewhere 



^ 2^, and refs. therein. The 



zeroth moment, < v'^{t) >, is related to the average number of molecular clusters in the system, and 
usually shows an algebraic decrease with time The first moment, < v^{t) >, is related to the total 
volume which is a constant value for closely-packed agglomerations 22] and an increasing function 



n 

of time for loosely-packed agglomerations 8], cf. Fig. 1. From the expressions of both moments, 
it follows that the average cluster radius, Ravit), behaves as a power law in time, with a growth 
exponent l/{d+ 1) that apparently contains some signature of random close-packing of matter by 
having included the super-dimension d + 1 2^. {d + 1 tells us something about the minimum 
number of non-overlapping neighbors of a given cluster in a d-dimensional space.) These constitute 
the main characteristics of the model agglomeration/aggregation process in its late-stage {t >> 1) 
limit. 

The question remains about asymptotic values of the moments < (t) > that must be known when 
app lying formula H17|l . For closely-packed agglomerations, the moments are found to obey a power law 
I22I 

<D"(i) >~i("-i)/(2-") (n = 0,l,2), (19) 

whereas for matter aggregation one finds another power law 

<'y"(i)>-tt("-i)+"l/(2-") (n = 0,l,2), t»l. (20) 

Notice, that for a = (c? = 1) both power laws above approach the same form, namely < i;"(t) > ~ 
i("-i)/2 When utilizing (|17() and (|19() it appears that for closely-packed agglomerations, cr^(t) can 
be fully identified with the inverse of < w°(i) > (the average number of clusters), cf. 2^ for details, 
what because of the constancy of < v^{t) >, leads to a^{t) oc Vsp{t), where Vsp = Vsp{t) ~< v^{t) > 
I < v'^{t) >, and can be termed the mean specific volume of the tightly-packed agglomerate, being 
equivalent to the inverse of its mean number density. The specific volume fluctuations read 

a^{t)oct^/^''+^\ (21) 

and if d — > 00, cr^(t) ~ Vsp oc t. 

When using H17|) and H2U|) . however, it turns out that for loosely-packed agglomerations cr^(t) is a 
quantity equivalent to the average cluster radius Rav (t) , see ^ lof . They behave in time as 

(T2(i) ocii/('*+i). (22) 

When d 00, a'^{t) ~ Rav{t) const, which means, that on average the system ceases to grow. Note 
that the standard diffusional regime, is always characterized by the one-half exponent, is achieved 



exclusively for the aggregations in c? = 1 because the only linear characteristic is Rav = Rav{t)- Note 
that Vsp{t) is not a linear characteristic, since Vsp(t) oc [Rav{t)]'^ usually holds. Here the d = 1-case 
must clearly be disqualified as standard diffusional, cf. Eq. 12f|) . 

Commenting on the last relations, 1)21(1 and (|22|) . one might furthermore conclude that they reflect 
a well-known Onsager conjecture that the fluctuations in a system undergo the same type of changes 
as the corresponding macroscopic dynamic variables : Here one may think of the specific volume of 
the agglomerate and the grain radius, and their behavior in the late-time domain, respectively. 

D. Coupling the instability (growing) and mechanical stress relaxation modes of matter 

agglomeration 

Poisson was likely the first who recognized that viscoelastic properties of fluids and solids can rea- 
sonably be compared in a suitable, mostly short-time domain, though the specification of the domain 
must be more precise for specifying the systems of interest. Maxwell successfully followed the ideas of 
his famous French predecessor, arriving at his well-known, in general non-Mar kovian, model of relax- 
ation 1^. In what follows we present our Maxwell-model-based ideas on how to distinguish between 
the two agglomerations under study, and how to switch on a kind of coupling between the (late-stage) 
growing and relaxational modes in the viscoelastic d-dimensional matrix that we investigate. The 
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and involves generically 



existence of the coupling seems to be experimentally justified, see 
the viscoelastic nature of the mega-cluster late-stage formation 

Thus, the afore presented rationale toward quantifying the fluctuations of the system can be strength- 
ened with a supporting phenomenological argumentation. The idea comes from a " coupled" diffusion- 
relaxation picture that appears in such a complex system. In any diffusion- migration growing process, 
the mechanical strain-stress fields play such a role as well. In our case, such a situation can be safely 
expected in the temperature domain T < Tpass ■ Another type of relaxation of the stress field, say am , 
is expected to prevail when the closely-packed agglomeration conditions are met. A different behavior 
may be observed when the closely-packed agglomeration conditions are lost for the first time, that is, 
at r = Tpass, when the loosely-packed context appears. In both temperature regimes, the relaxation of 
Cm (t) over the course of time, is very likely to go in a way essentially described by the current ^ . This 
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is expected to occur :22j presumably under (nearly) homogeneous strain conditions, e„j « const, for 
i >> 1. For an additional motivation of coupling matter agglomeration and stress relaxation picture, 
related to fracture phenomena, see ll 221 . 

n. 

From 22] it can be learned, that in the absence of non-Arrhcnius or fractal type kinetics, seemingly 
modifying the diffusion coefficient D{v) 2^, one expects the Maxwell dashpot-and-spring model to 
reflect properly the relaxation behavior. We wish to set up here a phenomenological picture, showing 
that both agglomeration and mechanical stress relaxation, where the stress relaxation takes place under 
slow growth conditions, proper of a mature growing stage in a viscoelastic multiphase medium 
are coupled processes |5|, . To work out the problem quantitatively, we will represent one of the two 
contiguous and matter-exchanging clusters in the agglomerate, say cluster (grain) 1, as an expanding 
one, equivalent to the spring, growing at the expense of its neighbor, to be named cluster (grain) 2, 
i.e the dashpot, to which, according to the Maxwell model, the contracting action should be assigned, 
c/., Fig. 2 for details; see Q. 

For the system with non-wide gaps, the Maxwell model conditions are almost satisfied, so that the 
two-cluster action can be extended over all pairs of contiguous clusters until the expanding (growing) 
eventually survive. In a next step, the same kind of competition appears as in the well argumented 
Laplace-Kelvin- Young scenario suitable for cellular systems [9| . This picture holds in the closely-packed 
context. 

In the loosely-packed context (a system with wide gaps), we may have qualitatively almost the same 

ii . . . . . 

picture |3 but with several differences which implies that cluster expansion would not be likely so 

vigorous. Since the corresponding gap is wider, therefore untight, the fluid leakage might be more 
pronounced. Thus, the fluid response against the piston wall is weaker, and the Maxwell type stress 
relaxation no longer applies, cf. the caption of Fig. 2. 

The stress relaxation can be described by introducing an exponent x in Maxwell-like quasi- 
fractional model presented here below. This exponent should be, in general, d-dependent, and points 

n 

to a difference when comparing with the classical Maxwell model 2:j. Here, we offer a coupled matter 
diffusion and stress relaxation picture, but for a random walk performed in the configurational space 
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FIG. 2: Maxwell sequential spring-and-dashpot (quasi-fractional) model with narrow (d-independent) and wide 
(d-dependent) gaps, shown schematically in two subsequent time instants ti and t2, where t2 > ti > 0, from 
left to right, respectively. Grain "1" consists of the spring and the piston's upper wall, to which the second end 
of the spring is attached, while its first end is mounted either on Oi or O2, from left to right, respectively. Grain 
"2" consists of the viscous medium inside the cylinder as well as the inner wall of the piston. The cylinder's 
walls complete the overall model structure of the viscoelastic grains. The material exchange between " 1" and 
"2" is assured by the existence of the gaps: narrow CPi_2 gaps in case of closely-packed agglomeration, and 
some two wider (here, represented by the left-hand side gap, LPi''*') in case of the aggregation. Therefore, the 
piston- and-cylinder system, containing a viscous fiuid, here composed of big and small particles, is either more 
(densely-packed agglomeration) or less (undensely packed formation) leakproof. The overall material exchange 
is caused by spring expansion along z axis, which results here in a longitudinal expansion of grain " 1" at the 
expense of grain "2", cf. \http://www.j-npcs.org/abstracts/vol2000no4-html\ Notice formally that: OiE = O2E, 
and for ti one has zi{ti) = Oi-E, 2:2(^1) = BE as well as for t2 one gets Zi{t2) — O2I, ^2(^2) ~ IE, which 
results in grain expansion-contraction behavior, like OiB < O2I and BE > IE, when mutually comparing the 
distances along the z axis at ti and t2, respectively 



2] 



22| . As is known, the Maxwell stress relaxation picture leads to an exponential decay of the stress: 



C^rnit) exp{~t/TM), 



(23) 



Einstein-Stokes-like formula 



9 



where tm is a reference time for the concentrated clusters ^ to be eventually inferred from the 

Q 2^ 



This behavior holds for T < Tpass- As mentioned above, for 



T > Tpass we propose 

d(7r, 



dt 



0, (24) 



where the above is usually true when the internal strain field, e,„ is practically constant, see above. 
When solving H24fl . one obtains 

o^{t)^{t/T)-^'^^-^\ (25) 

where x = 2(i + 3; about t, see Q or Q. Notice that for x = 1 and r = tm in Eq. one gets 
the solution H23|) ; for x 7^ 1 Eq. H25|) is the only solution to the relaxational problem as stated. When 
comparing Eqs. H23|l and (|25|l one sees that the relaxational response goes slower for the late-time 
loosely packed aggregational context than for its densely-packed agglomerational counterpart. 

III. QUALITATIVE SIGNATURES OF CHAOS IN MATTER-AGGLOMERATING 

SYSTEM 

Let us consider a few qualitative signatures of chaos in matter-agglomerating systems from the 
literature. Such certain signatures for systems of the type studied in the present work are summarized 
in Table 1. 

The items stated in Table 1 do not exclude other possible forms to chaos, or its signatures, in matter- 
agglomerating systems. We do not pretend to describe all of them, or even their majority. For routes 
to chaos recommended from physical point of view one would usefully consult j^^, 2^; which routes, or 
scenarios of chaos, are recommended by mathematicians, especially when a partial-differential-equation 

n 

formalism of F-P-K type is effective, can be found in and in refs. therein. 

IV. SOME QUANTITATIVE MEASURES OF CHAOS SIGNATURES IN 
MATTER-AGGLOMERATING SYSTEM 



In Q 



some entropic-like nonequilibrium measures of growth 



, (.) /ilVWl) , ,^1,2,3,... (26) 

^ Hi) ^ for t»l 



TABLE I: Qualitative signatures of chaos in a model matter-agglomerating system of interest 



Item 


Signature 


Refs. 


One: 


Entropic system of molecular-chaotic behavior 


[iiii, 24] 


Two: 


Lack of matter - depletion zones around the charged clusters 


[13, 22, 25] 


Three: 


Competition-and-loss effect: (un)tight spring-and-dashpot model 


[2, 23, 26] 


Four: 


IBCs of normal (e.g., Neumaim) or abnormal type should be taken appropriately [^^^, ^L] 


Five: 


Bethe-type frustration in coupled relaxation and late-growing events 


[L 2^, 2^ 


Six: 


Growth viz instability: random close-packing with its d + 1-account 


[a 2^, 24] 


Seven: 


Nonequilibrium entropy measures viz mean-harmonic speeds 


[2,i4, 28] 


Eight: 


Entropic potential(s) assuring nonequilibrium character of the phenomenon 




Nine: 


G-K type construction of D{v), B{v), and its consequences 


[2,i4, 27] 


Ten: 


Diffusion-space pre-chaotic (Fibonacci) feature by Da 


[S^ii, 22] 



as well as for the mechanical stress relaxation evolution 

(.) , rf. 1,2,3,... (27) 

^ ^"W ^ for t»l 

have been proposed. This seems to be working most appropriately in a growth-and-relaxation syn- 
chronization metastable regime 

^ Rav-'^^ (28) 

y 



which represents the H-P-G condition 



appropriate for the fluctuational late-time regime 



(29) 

3 



of interest here. Bear in mind that if certain 



empirical modifications of the formula (|29|l are applied toward obtaining a specific form, interconnecting 
Eq. H26() with Eq. (|27|l . one gets something like 



= qt^sp^^K (30) 



where typically q > 2. In the classical H-P-G limit J^sp*-'^-' — '^fisp^'^'' holds. However q may also strive 
for obtaining superplastic effects, i.e. when taking on fractional values, cf. 0|, and refs. therein. This 
is sometimes termed in physical-metallurgical literature the inverse H-P-G effect. 



Because of Eq. (|^ 



'-P^^^^dTT' '^^I'^'S,-- (31) 



Realize that formulae (|28() and 1)2911 might again be interpreted in terms of the Onsager conjecture 



see above. Since the overall exponent in 1)25(1 reads 

1 1 



(32) 



X-1 2(d+l)' 

which is exactly one half of the growth exponent Vgp'''^^ given in H22|) . see (|31|l too, one consequently 
provides 

^'^p^'^ = 2(dTT)' (33) 
Let us emphasize here that stands for the so-called Nutting exponent for relaxation, and can be 
interpreted in terms of the loss tangent, that means, a well-known dissipation factor in the relaxation 
phenomena, mostly in dielectric (e.g., macromolecular) environments (22l |. 

V. NUMBER-THEORETIC MEASURES OF SPATIAL AND TEMPORAL 
IRREGULARITIES IN AGGREGATION-AGGLOMERATING SYSTEMS 

It is interesting to note here that x — xi^) , i-S. 

X(d) = 2(rf+1) + 1. (34) 

A certain generator of the Bethe-lattice elements, is recovered starting from the 3-bond (initial) gen- 
erator for d = 0, and continuing with d, upon identifying d as the numbers of emerging bonds in a 



gelling system 



This is a very useful tool for the mean-field description of gels, and other multi 
bond-containing systems. T 

In this way, an odd number Bethe-lattice generator for subsequent d-s can be offered, see Fig. 3. 
Another, equally interesting observation can be offered, namely 

2 1 1 



1/ (2) 1/ (1) 1/ (3) ■ 
'^sp ^sp ^sp 



(35) 



This means that for loosely-packed agglomerations the harmonic-mean rule for the growth speed 
is exactly fulfilled. Mutatis mutandis, one can expect the same type of rule, Eq. ()35|) . for iisp^'^\ 




FIG. 3: An example of a not much developed Bethe lattice, which by itself manifests a frustration because of 

2^, 



" having problems" with containing all of its nodes in the available d-dimensional space \2(^ , to some extent so 
as, for example, the population of Verhulst fellow countrymen does in the available Belgium territory [2^ 

Let us recall that the fluctuations (T^(t) have been proposed as a reliable criterion of diff'erentiating 
between aggregation and agglomeration, and that an efficiency (harmonic-rule, see (|35|) ) additional 
criterion, derived from the constructed fiuctuational proposal, supports the aggregation in dimension 
d {d — 1, 2, 3), with an emphasis placed on c? = 2, where 'golden-ratio- like' or harmonic-mean properties 
are in favor. The mean-harmonic speed implies that the center of mass of a moving body, referred 
consequently to as the molecular cluster, may not span the same distance, say s, back and forth, during 
a time period. This leads to a quite realistic quantification of a mean speed on the distance 2s, and 
also shows that not an arithmetic mean of the back and forth speeds but a harmonic mean applies 
here. Such a schedule can likely be extended over the agg regation that essentially relies on random 
matter attachments and/or detachments of particles ^ 23, in such a way a forward sub- process may 
essentially go unidimensionally, like in a ballistic motion, whereas its reverse counterpart would explore 
the whole three-dimensional domain. This situation typically appears in the case of matter desorption, 
in which detachment occurs part by part from a "reactive" surface spot. 

Some other confirmation of (mean) harmonicity, its close relation to the golden rule, and to the 
Fibonacci sequencing (characterizing well the fractality of "diffusive" microstructures), is hidden in 
the (macro)ion (or, cluster [2^) diffusion coefficient in our model, Eq. Hl()|) . which is also included 
in the free energy <&. The label, or the lower index, is simply a, which for d — 1 results in a = = 0/1, 



for d = 2 gives a = 1/2, whereas for d = 3 offers a = 2/3. The values of the diffusion coefficient (Eq. 
H16|l 'l are: Dq = D/1, Di^2 = D/2 as weU as to 132/3 = -D/3, respectively. They correspond to the first 
five-number Fibonacci sequence, composed of the numerators and denominators of as, like 0,1,1,2,3, 
and obe a„+2 = Oji+i + Oji, for the three subsequent Fibonacci numbers a„, a„+i and a„+2- If so, one 
can provide the following two recursive formulae 

d^l, (36) 

and 

^(d-i)^^^ d = 2,3, (37) 
where ap = 0, ai = 1, 02 = 1, 03 = 2, 04 = 3 are the first five Fibonacci numbers. Since the analogy 



with gelling systems seems evident 



, this cannot be taken entirely as a surprise. The bonding in 



gels clearly goes as a branching process, being (as in the case of ultrametric space) quite naturally 
described geometrically in terms of Fibonacci numbers, thereby involving the notion of fractality ^. 

When finishing this section, let us note that both the characteristic chaotic measures, cf. Eq. (|31|l for 
example, have their random close-packing account d+1 involved. This is a landmark of randomness but 
readily appears as a space-filling action of modeled matter reorganisations. Realize that our rationale 
may apply just in the same vein to cluters-containing assemblies, evolving in a d-dimensional space, 
where a cluster is characterized by its fractal dimension 0<dF<d, cf. HQ 

VI. CHAOS IN AN INFINITE-DIMENSIONAL AGGLOMERATING AND/OR 

AGGREGATING SYSTEM 

Consider the case limc;_>oo. A corresponding chaotic measure for the late-stage growing event in the 
agglomeration of matter, very reminiscent of nonequilibrium correlational entropy measure reads 

= lim (M^) , (38) 

d^ooV ln{t) J 

whereas its counterpart for the relaxation is given by an analogous formula, namely that 

- lim (dl^) , (39) 



for t»l 



holds. They are consistent formaUy with the so-caUed correlational entropy (Kolmogorov-type) mea- 
sure, defined in 28] and follow the rationale presented in 3, in which some measures of chaos in 
dynamical systems described by partial differential equations have been discussed. For "thermostatic" 
systems out of equilibrium one has to speak of the so-called generalized fractal dimension formalism, 
first introduced by Grassberger and Proccacia, see [2^, and refs. therein. 

The most attractive reason for introducing such measures arises from the fact that if one evaluates 
both i^sp''°°-' and ij,sp^°°\ one unambigiously gets for the aggregation 

i^sp^""^ = - 0, (40) 

whereas for the close-packed agglomeration one provides 

'^sp^^^^ = 1, (41) 

and 

/i.p^°°^ = 00. (42) 

Thus, for both cases, Eq. H41|) and Eq. 142II . one arrives at a chaotic behavior in the nonequilibrium 
system 2, 2^ of a densely-packed agglomerate. This is not the case of the aggregation for which the 
common measure of its chaotic character is zero, cf. Eq. ()4U|I . 

Thus, proceeding consistently with the approach offered in we may conclude that the late-time 
aggregation process develops in an ordered way. The case d — 2 appears to be the most efficient since 
the harmonic-mean rule (|35|l is applied; for it the nonequilibrium character of the random process 

. n 

should be emphasized jSOJ. It is intriguing to realize that the system property called the harmonicity 
throughout is very much related to its nonequilibrium entropic or chaotic characteristic(s). 

VII. CONCLUDING ADDRESS 

Based on the above, we are allowed to state the following: 

(i) in matter-agglomerating systems chaos is revealed as a complex spatio-temporal and temperature- 
dependent phenomenon; 

(ii) nonequilibrium chaotic measures of any late-stage matter agglomeration modeled can be proposed 



relying upon the nonequilibrium Kolmogorov-type entropy measure, which makes a rehable (harmonic) 
quantification of the tempo of the process; 

(iii) coupUng late-stage matter agglomeration with relaxation of assisting elastic fields via an Onsager- 
type or, in the parlance of physical metallurgy, H-P-G conjecture ^, leads to several characteristic 

sub-effects (Bethe-lattice generator, first-five Fibonaci-number signatures, random close-packing d+ 1- 

. . J . . . . 

criterion |23l |. etc.) having their rationale in fundamental properties of the entropic or harmonic-mean 

character of the phenomenon 3^ ; 

(iv) as for the formal point of view: The presented mesoscopic system. Section 2, serving to describe 
the matter aggregation can be derived rigorously based on the Gibbs entropy production equation 



22| . and 

(v) its chaotic signatures can be inferred as presented in Sections 3-6, supported somehow by the 



ideas contained in [6|; at this point, a g eneral task remains to be done as to connect the type of chaos 



with the entropy-based scheme [10|,|2^ used to derive the equations of F-P-K 11], or diffusion, types 
^ ^ 27|, and how far the proposed measures of chaos (also, the ones used in the present review) 
are reminiscent of those used conventionally in nonlinear science ^J, [27, 2^ 3^? Perhaps, the 
Edwards' entropy measures for slowly moving grains, evolving (bio)polymer- or colloid- type matrices 
and compacted powders could also contribute to solve the problem |32l |. 
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